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I. INTRODUCTION 

The electron gas undergoes, in some conditions, a phase transition to bound pairs of electrons (the so-called 
Cooper pairs), that replace electrons as the fundamental agent of the electronic properties. The Cooper pairs are, 
from the point of view of their electric charge, objects qualitatively identical to the underlying electrons. From the 
point of view of their spin, on the other hand, they become integer-spin particles, that is, from the spin-statistics 
theorem, bosons rather than fermions. This shift of statistical paradigm of the carriers, from Fermi to Bose-statistics, 
results in the outstanding behaviour of superconductivity, that is, conduction of electric charge by a macroscopic 
coherent wavefunction (akin to a Bose-Einstein condensate). It has taken some time to capture the fundamental and 
universal features of this phenomenon and set them apart from particularities proper to certain cases only. The gap of 
excitations, responsible for zero resistivity, for instance, results from the long-range nature of the Coulomb interaction, 
but gapless superconductivity is also possible. One of the central, fundamental concepts of superconductivity is that 
of a coherent quantum state of charged bosons. Although superconductivity was discovered empirically, and its 
theoretical construction consisted in assembling a puzzle, it is now possible to envisage engineering superconducting 
phases in other systems, based on this understanding of condensation of charged bosons. If superconducting phases 
can be identified in other systems, progresses will be quick for the understanding of cuprate superconductivity, which 
still eludes compelling theoretical explanation of its intrinsic mechanism. 

A system that is making rapid and impressive progress in terms of creating and controlling macroscopic quantum 
states is that of microcavity exciton-polaritons [Ij (see [2] for a review) . These quasi-particles that combine properties 
of light (cavity photons) and matter (quantum well excitons) have been noted for their predisposition to accumulate 
in macroscopic number in a single or few quantum states |3J H] . They have many advantages from a practical point 
of view, such as their 2D geometry, which allows straightforward manipulation by lasers impinging at an angle, and 
their short lifetime, which allows continuous monitoring of the system, reconstructing its internal dynamics also 
by angle-resolved spectroscopy [5]. The pumping can be either coherent (driving states in parametric scattering 
configurations) [BHTO] or incoherent (with a constant flow of unrelated particles relaxing into the ground state ) [TTI - 
[TB] . In nitride systems, the formidable claim has been made of room temperature Bose-Einstein condensation pTlllSj . 
Recently, there has been great interest in propagation of polariton fluids [El [20] and their superfluid properties [21] , 
with reports of quantized vorticity [H] [23] and persistent currents [53] . 

These rising stars of macrosopic coherence have also been proposed to service another much sought after quantum 
phase at high temperature: superconductivity. Polariton condensates cannot conduct electric current themselves, 
being neutral particles. One of the proposed implementations involves "quatrons" (or quadrions) rather than po- 
laritons [25j . Quatrons are bound states of two electrons with a polariton. They remain bosons from spin-addition 
rules but carry an electric charge. A Bose condensate of quatrons would, through its superfluid propagation, exhibit 
superconductivity. To date, however, the existence of the quatron, predicted theoretically |25] , has not been found 
experimentally. Recently, we have approached the problem from another, and more conventional, angle, that of the 
BCS mechanism |26j with an important new feature: replacement of phonons by Bose-condensed exciton-polaritons 
in the role of a binding agent between electrons. In the present paper, we generalise the model proposed in [26) to 
describe a wide range of hybrid Bose-Fermi semiconductor systems where a Bose-Einstein condensate of neutral quasi- 
particles (excitons or exciton-polaritons) coexists with a Fermi sea of electrons. We show that indeed such systems 
are promising for observation of (high-temperature) superconductivity. Moreover, as the superconducting gap and 
critical temperature appear to be very sensitive to the concentration of bosons in the system and the latter may be 
controlled by direct optical excitation of exciton, light induced superconductivity in semiconductor heterostructures 
appears to be possible. In this work we closely follow the BCS approach, generalised and adapted to the case of 
superconductivity mediated by a Bose-Einstein condensate. Solving the gap equation in this case turns out to be a 
non-trivial problem, requiring careful analysis. 

A. BCS with a Bose-Einstein condensate as a binding agent 

BCS is a pillar of superconductivity theory, which relies on three main tenets: 

1. Instability of the Fermi sea, 

2. Existence of an attractive interaction, 

3. Condensation of charged bosons. 

These are the three insights that were mainly contributed by Cooper, Bardeen and Schrieffer, respectively, and that 
they could assemble into the BCS theoretical edifice and exploit to reproduce strikingly or predict successfully most 
of the superconductivity phenomenology. 
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The first point follows from Cooper's observation [57] that an arbitrary small attractive interaction between two 
electrons on top of the Fermi sea leads to a bound state (the Cooper pair) , thanks to the truncation of the momentum 
space for states with wavevector k > kp (Fermi wavevector). This is a general result, that follows from Bethe-Goldstone 
equation for the two-electron problem. 

The second point is the identification of an effective attraction between electrons that normally experience bare 
Coulomb repulsion. This attraction is attributed for conventional superconductors to an interaction through phonons, 
the Bardeen-Pines potential ;28J, that consists, vividly, of one electron wobbling the lattice at a first time, which affects 
another electron at a second time (and at a larger timescale, since the lattice dynamics is much slower than that of 
electrons) [29]. If the frequency of the lattice vibration is smaller than that of the propagating electron, the net effect 
results in an effective attraction [Leggett offers an insightful toy model of coupled oscillators to capture the essence 
of the interaction character (attractive or repulsive) [50]]. 

The last point is the so-called ECS state, which is a coherent superposition of paired bound states and which brings 
the two-particle Cooper effect to a collective behaviour of all electrons in the system [31j . 

With these three ingredients put together, the BCS theory is complete [3^. For our purposes, points 1 and 3 will 
be regarded as fundamental and well established features of (BCS type of) superconductivity. Point 2, that might 
appear a mere desiderata for point 1 to apply, leaves us room for identifying and designing new types of attractive 
potentials, optimizing the range of applicability and strength so as to obtain robust superconductivity (e.g., holding at 
high temperatures or high magnetic fields) or its manifestation in a new class of systems (in microcavities). Bardeen 
himself, with coworkers [33], investigated possibilities to engineer a more robust BCS in a bilayer structure where 
excitons replace phonons as mediators of the interaction. The idea of substituting phonons by excitons was pioneered 
by Little 34J and developed by Ginzburg [35' , who coined the term and theorized the possibility of high-temperature 
superconductivity, much before it came to fruition with cuprates 1361 . Those, however, exploit another (still unknown) 
mechanism different to BCS [37] . 

In the following, we revisit the Ginzburg mechanism, based on BCS, with emphasis on maximizing the strength of 
interaction between electrons, so as to maintain their binding, and therefore superconductivity, to higher temperatures. 
The text is organized as follows: in Section [llj we give a short overview of the BCS mechanism and the exciton 
(Ginzburg) mechanism, outlining the points of special interest in our case. In particular, we introduce the gap 
equation. In Section [111] we introduce our hybrid Bose- Fermi system configuration, its Hamiltonian and microscopic 
interactions, and the effective electron-electron Hamiltonian that results from a mean field approximation for the 
condensate and the usual Frolich transformation. We obtain the shape of the effective electron-electron interaction 
U, that wc find to be quite different in character to the Cooper (square well) potential. We also consider possible 
variations of our scheme, namely, a microcavity with a condensate of exciton-polaritons and a condensate of indirect 
excitons in coupled quantum wells [38] . The system of coupled quantum wells explored by several groups [39] [40] might 
be easier to realize and study (it does not need a cavity) and presents some interesting differences as compared to the 
polariton system. On the other hand, polaritons condense at much higher temperatures |18| . In Section IV, we study 
the gap equation for a Bose condensate mediated effective interaction. Because the potential is not positive-definite, 
the problem is not well-posed numerically. We propose a simplified potential and an approximate solution of the 
gap equation, that we motivate by studying its validity on well-established approximations. We obtain the critical 
temperature in this case. In Section jVj we give our conclusions and perspective on this new application of excitons 
and polaritons and discuss how to measure the effect experimentally. 



II. BCS AND GINZBURG MECHANISMS 



Superconductivity is a fundamental property of solids. At low enough temperatures, most metals superconduct. 
As the mechanism is rooted in quantum mechanics, temperatures are expected to be very low and indeed this is 
the case for all metals. After considerable theoretical efforts from various groups, a compelling theoretical model 
was assembled by Bardeen, Cooper and Schrieffer, the so-called BCS theory [53]. The model provides the critical 
temperature: 

fcBTc^ee-l (1) 

where, in conventional superconductivity, hcuj^ is the Debye energy, and g — Af(0)V with Af{0) the density of electrons 
at the Fermi energy and V the electron-phonon coupling strength. The Debye energy is, in good approximation, the 
maximum energy that can be carried by a phonon. The mechanism therefore relies heavily on phonons, as was 
realized empirically before the advent of BCS (good conductors, for instance, are bad superconductors, since V is 
small, or through the isotope effect, which correlates critical temperature with weight of the crystal atoms, and thus 
with resonance frequency). The exponential form and the presence of J\f{0) show that the effect is a collective one 
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involving all electrons, which have formed a new phase of matter that cannot be approached perturbatively by the 
independent electron pictures (since f{z) = e~^l^ has no Taylor expansion around zero). Based on the theory, and 
accumulated experience, it was widely accepted that critical temperatures would not exceed a few tens of Kelvins, 
since the Debye energy, which can be quite large in some systems (hundreds of Kelvins), is exponentially reduced. In 
all conventional superconductors A/'(0)V^ ^ 1 (the so-called weak-coupling regime). 

Ginzburg made the obvious but daring assumption that to achieve higher critical temperatures — crucial for technical 
applications which can easily be understood to be momentous — it is enough to find a system where O and/or g are 
increased. Replacing phonons by excitons, for instance, Ginzburg found that values 8 w 10'^, lO^'K as well as g ^ |, 
I are obtained, yielding temperatures of several hundred Kelvin [4lj . High critical temperatures have been later 
reported in cuprates [36] and nowadays, temperatures as high as 125K are reported in systems such as Tl-Ba-Cu- 
oxide. Just as superconductivity in metals was discovered ahead of theory, cuprate superconductivity does not appear 
to follow the ECS pattern [32] even with substitution of the mediating field as proposed by Ginsburg. All attempts 
to date to realize the exciton mechanism or a variation of it have remained fruitless, but it is this mechanism, rooted 
in BCS, which we consider in this text. Before we turn to the mechanism we propose, which consists in substituting 
the phonons of conventional BCS, or the excitons of Ginzburg scheme, by a Bose-Einstein condensate (BEG) of 
excitons or exciton-polaritons, we recall that the starting point of a microscopic derivation of superconductivity is the 
gap-equation: 
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where C/kk' is the effective interaction between electrons with wavevectors k and k' and energy E. The gap A can be 
identified as the macroscopic wavefunction of a Gooper pair, which is also the order parameter for the superconducting 
phase. If it is nonzero, the system is in the superconducting state. A realistic microscopic treatment of U is very 
complicated. A simplified version is provided by the Jellium model, which is a toy model of a metal that gives 
predominance to electron-electron interactions, that is, in particular, the underlying crystal is approximated as a 
uniform (structureless) background (like a "jelly") in which the interacting electron gas evolves under its own self- 
interactions and the overall charge cancellation of the background. A popular effective interaction is given by the 
Bardeen-Pines potential |30j which is derived from the microscopic form of the electron-lattice interaction: 



(3a) 



(3b) 



(1 + ;^)(1-^) 



with q = k — k'. In Eq. (3a I, q-jp is the Thomas- Fermi screening parameter and Wph is the phonon dispersion. The 
first term between the brackets is bare Goulomb repulsion, and the second term, which is frequency dependent, follows 
from the perturbative coupling to the lattice, tracing out the phonons. 



The physical meaning of Eq. ( 3a I is at the heart of the phonon-mediated mechanism. The interaction is oj dependent, 
which means, in Fourier transform, time dependent. This reflects the famous retardation effect in superconductivity. 
This effect is based on the strong difference between the electron Fermi velocity in metals and the sound velocity. 
Roughly speaking, a fast electron from the Fermi surface creates a slow phonon and goes away. After some time, 
another fast electron arrives and absorbs the slow phonon. The average distance between these two electrons remains 
of the order of lOOnm, the distance at which the screened Goulomb repulsion can be safely neglected. Due to the 
retardation effect, a weak phonon-mediated attraction of electrons wins over their Goulomb repulsion and provides 
formation of Gooper pairs at low enough temperatures. 



III. INTERACTION HAMILTONIAN 



A sketch of the structure we propose appears in Fig. [T|^c) . A QW is doped negatively (with density of electrons n) 
and is put in contact with another QW where excitons are formed and are stable. The first (lower) QW hosts the 
2DEG which is to undergo superconductivity while the second (upper) QW hosts the condensate of excitons that is 
to mediate it. This structure can also be placed at the maximum of the optical field confined in a microcavity |26| 
(formed by two Bragg mirrors facing each other) . The latter configuration also has the advantage that two excitonic- 
QWs can sandwich the 2DEG-QW since in this case the condensate is delocalised in the entire structure thanks to the 
photon fraction of the polariton. This allows an increase by a factor 2, or more if the multi-layer structure is further 
repeated, of the densities achievable in this system. The drawback of microcavities is a short radiative lifetime of 
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FIG. 1. One of the possible designs to evidence exciton-polariton mediated superconductivity, (a) Polariton (thick blue solid 
line) and exciton (purple dashed line) dispersions with schematic representation of scattering of the boson that mediates the 
interaction between electrons. Also shown in thin blue solid line is the renormalised dispersion -Ebog, essentially identical to 
Epoi- (b) Corresponding scattering of two electrons on the surface of the Fermi sea, exchanging momentum q through scattering 
of a boson in panel (a), (c) detail of the sandwich structure, showing the doped well (below) containing free electrons and the 
well hosting the condensate (above) of excitons or of exciton-polaritons. This layer can be sandwiched in different ways. 



polaritons. In order to maintain the polariton condensate one needs to pump it resonantly by a high intensity laser, 
which may lead to undesirable heating of the system. The alternative system which we consider is a condensate of 
long-living spatially indirect excitons in coupled QWs, separated by a thin but high barrier from a QW containing a 
two-dimensional electron gas (2DEG). This system is simpler to realise as it does not contain the optical cavity. The 
condensate of indirect excitons may be maintained by a low intensity optical excitation. Moreover, indirect excitons 
in coupled QWs possess very significant dipole moments which strengthens their interaction with electrons from a 
2DEG. However, Bose condensates of indirect excitons have been found only at temperatures below IK [38j, while 
condensates of exciton-polaritons are now routinely produced at room temperature |43j . In the following, we compare 
the advantages of the two respective schemes. Most of the underlying model applies to both equally. We focus more 
particularly on the microcavity system since it is more general. The exciton system is recovered by decoupling the 
photons. 

In Fig. [ija), the exciton (resp. polariton) dispersion is shown in dashed purple (resp. solid blue). The difference 
between the bare polariton dispersion (thick blue) and the bogolon dispersion (thin blue) is very small over the range 
of exchanged momenta of interest (of the order of the Fermi wavevector). The condensate in both cases is at fe = 0. 
Scattered particles at any wavevector between — 2fcF and 2fcF mediate electron-electron interactions on the Fermi sea, 
as sketched on Fig.lljb). The model microscopic Hamiltonian is taken as: [26] 



H = 



|^Spoi(k)a|^ak + Eci{'k)al(7k 



E 

ki,k2,q 



Vc(q)'7ki+q'^k2-q'^kiO-k2, +^^x(q)crki^ki+qak2+q"k2 + ^L+qOki-l-qaka 



(4) 



with £'poi(k) and i?ei(k) the polariton and 2DEG dispersions for the in-plane wavevector k, respectively. In the 
exciton case (without the microcavity), it suffices to replace Epo\ with Ecx in the above. Vx is the electron-polariton 
interaction, U the polariton-polariton interaction and Vq the electron-electron repulsion. We now consider these terms 
in turn. Parameters assumed are listed in table HI 



A. Electron-electron interaction 



In the original BCS mechanism, electron- electron repulsion is either neglected altogether or overcome by the at- 
tractive mechanism and not manifested outside of the attractive window. We take it into account here since it is a 
detrimental factor for binding and most of our concern for experimental realization is to optimize this value. The full 
form of the potential is given by the Yukawa potential: 



(5) 
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Value 


e 


Permittivity 


7eo ~ 0.2 As/(mV) 


Pe 


Electron reduced mass 


0.22 ^ 1 ic 
0.22+1.25 




Hole reduced mass 


1 25 n or 

0.22+1.25 


L 


Distance between wells 


5 nm 




Coulomb screening length 


^i.z X iU m 


mx 


Exciton mass 


(V.2Z + i.2ojme ~ i.o X iU Kg 
iU me ~ y.i X iU Kg 


mc 


Photon mass 




R,abi splitting 


A ^ Tvip\/' rnl 
40 me V [uj 






1 /a/2 fll 
1/ V z [J-J 


ftp 


Fermi wavevector 


5 X 10* m~^ 


OB 


Exciton Bohr radius 


1.98 X 10"^ m 


Ry 


Exciton Rydberg 


32 meV 


I 


Dipole moment 


4 nm 



TABLE I. Parameters used in the numerical simulations. In square brackets, the value for the exciton case when the parameter 
differs from the polariton case, otherwise parameters have been taken the same for comparison. 




3D 2D 

FIG. 2. Average over the Fermi Surface, in 3D and 2D. 



with screening constant k. We get rid of the momentum dependence by averaging the potential 14ff(aj,q) over the 
Fermi surface (FS), where: 

q = ki - ka (6) 

and ki 2 are the momenta of the two interacting electrons on the FS (so that, in particular, ki = k2 = fcp), as shown 
in Fig. jl] In 3D, the FS is the surface of a sphere, while in 2D, it is a circle (we shall speak of surface in both cases). 
The vector difference is therefore joining the two end-points on the surface. If the potential has spherical symmetry 
{V{q) = V{q)), the average of all two vectors on the FS reduces to that where q is pinned at one point of the surface 
(the south pole in Fig. [2]) and runs overs the FS. This average, in the particular choice of Fig. [2| is the usual polar 
integration with k2 describing the surface as 6 (and (j) in 3D) are varied, with: 

q2 = 2fc|(l + cos6'), (7) 

from Al Kashi's theorem, so that the average potential Vee reads, in 3D: 



/'27r PTT 

Zsiuj)^ / Voff{q,uj)kl sine dOd^/J^, 

JO JO 



(8) 
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2k-p/ 

FIG. 3. Average electron-electron repulsion (in natural units) in the QW. This should be made as small as possible, which is 
obtained for higher Fermi wavevectors (for a given screening length). 



where M is the normalization, i.e., the same integral where is replaced by unity. This gives, in 3D: 



T4ff(^/2fc|(l + cos6'),a;)dcos( 



t/cff(fcF%A?,w)dl?, 



where we integrate over d — 2(1 + cos 6') since this is a natmal variable in Eq. ([t]), and, in 2D: 



1 

2^ 



27T 



V^s{j2kl{l + cose),uj)d6 . 



Our system is 2D, in which case, from Eqs. ^ and (11): 

2 TT — 2 arctan 



Vc 



2tt 



2fcF 



d9 



47reA Jo v/2/s|(l + cos 61) + k 



e 

2tt€A 



In 



g2 l-2feF ^4fe|-K2 



if K > 2/cfi 



if K < 2A:f- 



(9) 
(10) 

(11) 



(12) 



(13) 



This is plotted in Fig. [3| Since we are trying to maximize attraction, that is, minimize repulsion, systems with 
small screening length and large wavevectors should be favored (but these parameters play critically on other aspects 
of the mechanism and the optimum is not compulsorily kp/n 3> 1). 



B. Electron-exciton interaction 



The electron-exciton or exciton-polariton interaction is one of the most important ingredients of the mechanism, 
as it ultimately determines the shape of the effective potential. In the microcavity, an electron (from the 2DEG) 
interacts with a polariton (from the condensate) through its excitonic component, so this is really the electron-exciton 
interaction that is to be computed, weighted by the Hopfield coefficient (the excitonic fraction) X. Let us consider, 
therefore, the scattering of an electron in one of the parallel QW, separated by a distance L from the QW with 
excitons. The matrix element of the direct interaction between excitons and electrons reads: 



Vx{(l) = / n(Q,re,r,0«'*(k,ri)V^(ri,re,r/,)^'5,(Q,re,r^)**(k,ri)dridredr,, 



(14) 
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FIG. 4. Electron-exciton interaction in the geometry of Fig.jTJ decomposed as the direct interaction (dashed purple) and dipolar 
interaction (solid blue) when the exciton is induced with a dipole moment d. The latter is both much larger and maximum at 
zero exchanged momentum. 



where ri,re,r/i correspond to the 2D coordinates of the 2DEG electron, the exciton electron and the exciton hole 
respectively. The 2DEG electron is described by a plane wave while the electron/hole in the condensate are assumed 
to be in the Is bound state with plane wave center-of-mass motion: 




mQ,re,r,,) = W A J_t7^(^^)[/,(^,Je*Q-(/3er.+/3.r.)^-|r.-r.|/a. ^ / 2 C^feH^M ^^Q • Rx /a. , (Igb) 

V ttA qb V ttA aB 

where Rx = l3e^e + PhJ^h,^x = Te — r/i are in-plane coordinates of the center of mass of the exciton and relative 
coordinate of electron and hole in the exciton, Ue (ze) and Uh (zh) are normal to the QW plane electron and hole 
envelope functions, respectively. We also consider the existence of a dipole moment d for the exciton, which can be 
intrinsic to the structure, because of spatial separation of electrons and holes in coupled QWs, or (in the case of 
microcavities) be induced by an internal piezo-electric field, or result from an externally applied electric field. To 
account for all these possibilities, one can consider the layers of electrons and holes in the exciton shifted in the 
z-direction with respect to the position of the center of mass by a distance I — d/e L. The matrix element of the 
interaction is then computed to be: 



Vdiriq) 



-qL 



'2eA q 1 [1 + (/3,gaB/2)2]3/2 [I + {/3^,qaB/2rf' 



2eA 



[1 + {l3,qaBl2Yf^ [1 + {PhqaB/2Yf'^ , 



(16a) 
(16b) 



where Eq. (16a) is the direct electron-exciton interaction that exists even in the absence of a dipole moment of the 



exciton, and Eq. (16b) is the dipolar interaction. The direct interaction vanishes at small exchanged momenta, while 



the dipolar- induced one assumes its maximum value here of 2d/(2eoeA). Overall, the dipolar interaction is naturally 
much larger than the direct one, since the exciton is electrically neutralj^ These facts are summarized in Fig. |4] 



^ Note at this point that there is an error in Ref 1261 where only the direct exciton interaction has been taken into account with an 
incorrect power in the parenthesis, which led to an expression similar to the dipolar interaction. The direct interaction by itself turns 
out to be too small to evidence superconductivity with the parameters chosen in Ref. |26| . therefore a dipole moment should be induced 
in this case, say by applying an external electric field, to restore the effect. 
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C. Polariton-polariton interaction 



We treat the polariton-polariton interaction within the s-wave scattering approximation with strength U ~ 
Qa^RyX'^/A (where Ob is the exciton Bohr radius, Ry the exciton binding energy and A the normalization area. 
X is the exciton Hopfield coefBcient, the square of which quantifies the exciton fraction in the exciton-polariton 
condensate) [33]. For interaction between bare excitons, X — 1. Exciton-exciton (polariton-polariton) interactions 
are repulsive, in general. They result in linearization of the elementary excitation spectra, the Bogoliubov dis- 
persion, but its role is not crucial to the mechanism of superconductivity we discuss, since at the wavevectors of 
interest, the changes brought by this term are very small compared to the kinetic energy of non-interacting excitons 
(exciton-polaritons) . 



D. Effective interaction 



We now proceed to bring our microscopic model towards a form suited to study Cooper-pairing and superconductiv- 
ity, that is, we apply the canonical Frohlich transformation that will result in an effective ECS hamiltonian. Just as in 
the case of phonons, we start by getting rid of polaritons. We assume a condensate is formed with mean population Nq. 
We do not consider which mechanism, coherent or incoherent, is responsible for creating and maintaining this state. 
We do assume, however, it is coherent and with a definite phase, so that we can apply the mean-field approximation 
°'ki+q'^ki ~ ('^ki+q)'^ki + '^ki+q('^ki) and (ak) ~ \/NqA5\^_o with TVq the density of polaritons in the condensate. This 
allows us to obtain the following expression for the Hamiltonian, after diagonalizing the polariton part by means of a 
Bogoliubov transformation (that leaves the free propagation of electrons and their direct interaction, Hq, invariant): 

H = Y, £;ci(k)4ak + EhosO^)bib^ + Hc + J2 *^(q)4^k+q(fe^-q + &q) (17) 

k k k,q 

where i?bog(k) describes the dispersion of the elementary excitations (bogolons) of the interacting Bose gas, which is 
very close to a parabolic exciton dispersion at large k: 



£;bog(k) = ^i;poi(k)(£;poi(k) + 2Un^a) , (is) 

where £'poi(k) = i?poi(k) — i?poi(0) and with the renormalized bogolon-electron interaction strength: 

M(q) = ^N^XV^{<i)f- ' . (19) 

W V ^^^^Y2C/A^o^-ii;bog(q) + i?pol(q) ^ ' 



The last term of Eq. (17) coincides with the Frohlich electron-phonon interaction Hamiltonian, which allows us 
to write an effective Hamiltonian for the bogolon-mediated electron-electron interaction. This results in an effective 
interaction between electrons, of the type X^ki ka q ^ff('i' '^)''^ki''^ki+qfk2+qO'k2 ■ The effective interaction strength 
reads V;ff(q, w) Vc(q) + VA(q,a;), with: 

2Af(q)^i;;bog(q) 

l/A(q,w) = 1—-^ ■ {20) 



Equation ( 20 ) recovers the boson-mediated interaction potential obtained for a Bose- Fermi mixture of cold atomic 
gases [iSj, in the limit of vanishing exchanged wavevectors. It describes the BEC induced attraction between electrons. 
Remarkably, it increases linearly with the condensate density iVo. This represents an important advantage of this 
mechanism of superconductivity with respect to the earlier proposals of exciton-mediated superconductivity [331 
1341 146] , as the strength of Cooper coupling can be directly controlled by optical pumping of the exciton-polariton 
condensate. The attractive potential is displayed for various exchanged energies in Fig. [Sj as a function of the 
exchanged momentum expressed directly through the angle 6 defined in the Fermi circle (cf. Fig. [2]) . As commented 
earlier, the negative part corresponds to attraction, and the potential alternates between repulsive and attractive 
character, obtained at different exchanged momenta. We do not want to keep track of such complicated wavevector 
dependence, and therefore will average the interaction over the Fermi sea. A notable feature is, for most values of w, 
the presence of a pole 9o, where i?bog(\/2fc|(l + cos 6*0)) = hoj. As is seen in the figure, separates the attractive 
part from the repulsive part. The average will bring the additional convenience of cancelling such divergencies. We 
note as well that the spectrum of excitations of the exciton BEC may be changed in the presence of the electron gas. 
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FIG. 5. Effective electron-eiectron interaction as a function of tlie exchanged momentum q — v^'sf Vl + cos 9 on tlie Fermi 
sea at energies hoj = 40, 50 and 60meV, respectively. Tlie potential is symmetric around n. In the first case, the whole 
dependency (over [0,27r]) is shown, with a first zoom in the inset showing two poles and another the attractive region that 
is of small amplitude but extends over a large range. With increasing energies, the poles recedes towards smaller values of 9 
with a dominating effect of the repulsive (positive) energy. In the central panel, in dashed purple, the regularized potential is 
superimposed. 



so that their eventual dispersion may be different ii^bog [57". This has no effect on the Cooper pairing of electrons 
which we discuss here. 

We therefore wish to perform the average 

Va= / VA{q,uj)de, (21) 



where q — y/2kp(l + cos 9), as seen previously. Since Va is symmetric around tt we perform the integral only. 
The integral would be easily computed numerically if there were no pole. There are dedicated numerical methods to 
compute principal values numerically ^48, ,49^ , but in our case, since the pole is first order, it is enough to isolate it 
analytically by defining 

/(0) = (0-0o)^A(q(e),w), (22) 

so that then 

Jo t' — fo fo 

where /(^o) — linie-i-eo f{9) and the first integral is regular (it is shown in Fig. [s] in dashed purple). The integration 
(average) is then straightforward and produces the results shown in Fig. (|6]). 

If we take instead of i?poi a quadratic dispersion for the excitation 

E^k) = ^ (24) 

and assume all-excitonic interactions, X = 1, we can consider the same effect in the absence of a microcavity, relying 
on a purely excitonic (rather than polaritonic) BEC. In this case, the same procedure as detailed above leads to an 
effective potential as shown in Fig. [7] We kept all parameters the same for comparison except for the dipole moment, 
which we have taken three times as large {d = 12nm). This corresponds to the spatial separation of electrons and 
holes in the system of indirect excitons studied by Butov et al |38j. The potential we obtain in the exciton case is 
very different in character from the polariton case, and is more closely related to the Cooper (conventional) shape of 
a square well, or the Bogoliubov potential including repulsive elbows. 

From these potentials, one can proceed to solve the gap equation. 
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FIG. 6. (Colour online) Exchanged-momentum averaged interaction Uo{ijj) between electrons of the 2DEG, as A^o is increased. 
It is dimensionless and is attractive (resp. repulsive) when negative (resp. positive). Densities A'^o are shown in units of 10^^/cm^. 
In insets (a) and (b), a zoom of the regions delineated on the central (left) panel, firstly in the small energies (long time) range, 
where the attraction is seen to be repulsive for smaller densities, because of the Coulomb repulsion (shown in red), and secondly 
for the densities ~ lO'^^/cm'^ in the area where the character of the interaction changes abruptly from attractive to repulsive 
(higher densities are shown in dotted lines). 
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FIG. 7. Same as Fig. |6] (also in units of 10^^/cm^) but for the case of an exciton condensate. A larger dipole moment d has 
also been assumed. The potential is different in character, much closer to Cooper's potential with larger attraction at longer 
times. 
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FIG. 8. Numerical solution to the gap equation with the BCS square well potential (a). The gap A(^) is not stepwise 
as approximated in the BCS model, but its value at zero exchanged energy, A(0), is in close agreement with the analytical 
expressions, as seen in (b) and (c): A(0) as obtained from numerically solving the gap equation (dashed blue) and from Eq. (271 
(solid magenta), (b) in the weak-coupling limit when Uo <^ 1, with small deviations as coupling is increased, and (c) in the 
strong-coupling where the agreement becomes perfect again. Note that A increases linearly with the coupling strength out of 
the weak-coupling limit. 
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IV. GAP EQUATION 
A. Cooper potential 

The BCS gap equation ^ is easier to tackle as a continuous equation: 

J-oo 

where we have also introduced a finite temperature T from the Fermi-Dirac distribution of elementary excitations I50j . 
With the BCS approximation of a step potential, the gap equation at zero temperature simplifies to A(^,T) = 
— C/q J^^I^^^ ^{C'iT)/2E dS,' ■ If A ^ fi-^D, the ^ dependence in the integral boundaries can be neglected (or, coming 
back to Eq. ([2]), one sees that in the initial gap equation, Ak is exactly constant if ?7kk' = Uq and is zero otherwise). 
Thus, we can assume the gap to be of the form: 

^(^)^fA(o) ifiei</.-D, 

^' 10 otherwise. ^ ^ 



In this case, simplifying A(0) on both side of Eq. (25 1, we obtain 



A(0) = -fiwD/sinh(-l/[/o). (27) 



Equation (27 1 is better known as its approximation when [/q ^ 1, in which case it takes the form of the famous 



BCS gap expression, A(0) = 2fia;D exp(— l/i7o). 

Solving exactly the gap equation calls for some numerical method. Equation (251 is a nonlinear integral equation, 
of the type studied by Hammerstein, i.e., 

A(e,T)= / i^(e,e')/[C',A(0]dC' (28) 



where in our case if (x, y) = —Uo{x — y)/2 and /[y, z] — z ta.iih{^/y^~+^) / ^/y^ ^ z^. There are strong conditions of 
existence of nontrivial (nonzero) solutions when J/q < [51. 52^ , however the case when the kernel K is not positive 
definite, that is, in presence of repulsionj^ has been much less studied mathematically [S3]. Cooper's potential, 
being always negative, falls in the category of potentials which admit a unique nontrivial solution, with a stable 
numerical technique to obtain it, namely, since the mapping is contractive, by iterations of the gap equation: an 



initial (nonzero) function Aq is used to compute the rhs of Eq. (25), providing Ai which is injected back until the 



function converges. The gap computed in this way for a given Cooper potential is shown in Fig. 8[a). As one can 



see, the BCS approximation (261 remains a rather coarse approximation, since the gap turns out in this case to 
be bell-shaped rather than being a step function. Surprisingly, A(0) is however in much closer agreement with its 
approximation, as shown on Fig. |8jb) and (c) in the weak and strong coupling regime, respectively. There are small 
quantitative deviations in (b) between the numerical points and the formula when using the exact parameters of the 
potential. By fitting the numerical results, a perfect agreement can be found for slight variations of /iwd and Vq. The 
BCS approximation therefore turns out to be an exceedingly good one as compared to an exact solution of the gap 
equation. For the procedure to make sense, similar results should be obtained for a smoothed well that approximates 
the BCS square well [HJ. We will not address this point here, but go directly to the case where the potential is not 
always attractive. 



B. Bogoliubov potential 

The potential is not always attractive when, for instance, some overall repulsion, such as direct Coulomb interaction, 
is superimposed on the attractive Cooper potential, as shown on Fig.[9j The Coulomb interaction is time independent 
and should therefore extend to all w but here also a cutoff loq is introduced to avoid divergencies. This results in 
an attractive. Cooper-like potential, fianked by two repulsive windows. Such a potential is known as the Bogoliubov 
potential [511155] . 



^ as well as attraction, since the case of only repulsion admits only A = as a solution. 
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FIG. 9. Adding an overall repulsive Coulomb repulsion Vc (red, left) until a cutoff cjc to the BCS potential Vb (green, left) 
results in the Bogoliubov stepwise potential (right). We take the convention Vb, Vc, Uo positive in the above representation 
and Uo ~ Vo — Vc, so that, e.g., Vc < means attractive contribution of the "Coulomb repulsion". 



This approximation has been used to show extremely counter- intuitive behaviour of the gap equation and justify 
a-posteriori another heavily criticized approximation of BCS, neglecting Coulomb repulsion: the BCS mechanism 
indeed assumes only attraction between electrons, which can be dimmed by Coulomb repulsion, but which never 
explicitly appears as such (like in the Bogoliubov scenario). The great result of Cooper was that binding occurs at 
arbitrarily small attraction. An important result of the Bogoliubov potential is to show that the detrimental effects 
of Coulomb repulsion are greatly reduced in the gap [55] . 

We now present a linearization of the gap equation that allows one to obtain an approximate solution for the critical 
temperature [54 . By assuming the gap equation to be a two step valued function A = (Ai 
becomes (/ — 1)A — with 



A2)"^, the gap equation 



I = - 



Veil Vcl2 



where 



° tanh(e/(2fcBTc)) 



'"'"^ tanh(e/(2fcBrc)) 



, « In 



1.13/kt;]3 



d£, w In 



ujc_ 



(29) 

(30a) 
(30b) 



appear invariably in the columns of /. In Eq. (30a), the BCS approximation has been applied while in Eq. (30b I, 



the fact that |^| ^ has been used to neglect A in the denominator and the temperature in the numerator. The 
parametrization of the matrix in terms of the potential depends on the particular configuration (for example the relative 
widths of the various layers of the structure) . Here we have adopted the original parametrization of Bogoliubov, which 
assumes narrow repulsive elbows surrounding a large attractive central region. Solving the linear equation for Ii and 
then for the critical temperature, we find: 



ksTc ~ 1.13/iWDexp 



V 



1 



1 



Vc 



(31) 



In this form, one can see how Coulomb repulsion indeed introduces a small correction to the original BCS formula. 
The expression also seems to indicate that could be repulsive (Vb — Vc < 0) and still lead to a gap as long as the 
denominator in Eq. (31) remains positive. 



In the following, we compare these predictions with numerical solutions of the gap equation, keeping in mind that 
the iterative procedure is not assured, mathematically, to converge. We have observed that indeed, it sometimes 
encounters problems and exhibit strong instabilities, with bifurcations of solutions, for example. 

we show the evolution of the gap function as Vc is increased, from zero (BCS) to a point where the 



In Fig. 10 



overall potential is essentially repulsive. With onset of the repulsion, the gap acquires two negative sides and becomes 
a highly distorted function for large Vc . Note that the approximation of constant gap over the various regions is at 
least as good as for the case of BCS. 

In Fig. [TT] we now show the case where U q — — Vq is held constant as the strength of the repulsion is varied 
independently. In foresight of what is to come later, we also allow the elbows to be negative, that is, to contribute 
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FIG. 10. Gap (thick magenta and, magnified, thick khaki) of the Bogohubov potential (filled blue) solved numerically for the 
parameters: huu = 1, hujc — 2, Vo ~ 1 and Vc taking values from (a) to (d) of: (BCS), 0.1, 0.3 and 0.495. Coulomb repulsion 
results in a dip in the gap function that, for increasing values, results in oscillations in the gap function. Paradoxically, even 
when it is large and dominating attraction, repulsion does not prevent a gap, as seen in (d). 



an additional attraction to the conventional mechanism (for now we do not consider physical justification of this). 
Another unexpected result is obtained: the repulsive potential is in this case favouring a larger gap, as can be seen 
by comparing (a), where the gap function is highly oscillatory, to (d) where it recovers the BCS bell-shape. In (h), 
the distorted but overall attractive potential still results in a BCS type gap, but wider and larger. Note how the 
repulsion, by "squeezing" the gap, allows it to achieve much higher values than for the case of smaller or no repulsion. 

These unexpected results are confirmed phenomenologically by the Bogoliubov approximation Eq. (31), which 
we plot as a dashed line in Fig. |11[ compared to A(0) computed numerically. Here we should emphasize that the 
two quantities are not meant to be compared quantitatively (as we did when comparing the BCS formula with the 
numerical solution), since one, A(0) (computed numerically) is the gap at zero temperature while the other, fcBTc is 
the temperature at which the gap vanishes. There is a monotonous relationship between the two, that is, increasing 
A(0) implies increasing Tq, so one can appreciate the consistency of the results by observing similar trends. The 
obstacle to conducting an extensive numerical comparison is that it is an intensive task numerically to compute Tq, 
since this requires solution of the gap equation for various temperatures until the curve A(T) is obtained and its 
intersect with zero is found. A critical slowing down phenomenon makes the iterative process slower as the critical 
temperature is approached. In addition, numerical instabilities are stronger at nonzero T. Therefore, although it 
is relatively straightforward to compute A(0) numerically, it is not convenient to use this method to obtain Tq- 
On the other hand, the Bogoliubov approximation gives a fair estimate of Tq, but is not able to provide the gap 
at zero temperature, since at the core of its method, there is an assumption of vanishing A. Therefore, we have 
two complementary methods, each suited to provide a relevant aspect of the problem. We note that the gap at 
zero temperature is an important quantity which can be measured independently from Tq by Andreev reflection in 
conductivity experiments, which is why we discuss it here in a great detail. 

We can see, indeed, that the qualitative agreement is good and that the counter-intuitive features of the gap 
equation are reproduced by the analytical formula. 



C. Polar iton potential 



In the case of the polariton problem, we have seen that, even when neglecting Coulomb repulsion (as in the original 
BCS formulation), the potential U{uj) departs strongly from the Cooper potential and features two large attractive 
regions far from small energies, immediately followed by two strong repulsive windows. We extend the Bogoliubov 



method to a three-step approximation of this potential, such as displayed in Fig. 13 with, in reference to previous 
potentials, notations wd, wc and lub for the boundaries of the central, shallow attractive region, narrow, deep attractive 
region and repulsive region, respectively. This is a notation only and is not mean to be understood as referring to 
Debye, Coulomb or Bogoliubov in any strict sense. Following the same premises, we approximate the gap equation 
by a three-step valued function A = (Ai, A2, A3)-'". This approximation turns out to be an exceedingly good one in 
certain cases, such as the one displayed in Fig. |13[ Here there is even more room to choose a parametrization of the 
/ matrix. We now give general guidelines on how to build this matrix. The simplest method is to fix ^ on the Ihs of 
Eq. (25) at the center of each region and, in the corresponding row, take for each column the potential that is sampled 
more by the difference |^ — ^'|. Refinements are possible, such as weighting elements of / by coefficients which refiect 
how much time the variables ^ and ^' spend in the regions that determine the matrix equation. This problem has 
the following mathematical expression: how is the random variable X — Y distributed when X (resp. Y) is uniformly 
distributed in an interval [gi,gi^i\ (resp. [/i^, /ii_(-i]). The solution is easily obtained as proportional (normalize to 
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FIG. 11. Gap (thick magenta and, magnified, thiclc kliaki) of the Bogoliubov potential (filled blue) solved numerically for the 
following parameters: fioju = 0.5, huJc = 2, Vo ~Vc ~ 0.5 (fixed) and Vc (the height of the elbow) taking values from (a) to (h) 
of: -5, -0.5, -0.1, (BCS), 0.3, 0.5 (BCS), 1.5 and 5. The repulsive nature of the elbows change the character of the gap from 
bell-shaped to an oscillating function. The oscillating gap in the presence of strong repulsion allows, paradoxically, high values 
of the gap at zero-exchanged energy. 
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FIG. 12. A(0) for the case of Fig. [To] (a) and |11| (c), for their respective parameters, and (b), with hcj-o = .5, Vo = 1, uJc 
changing as indicated and Vc changing such that the area of the repulsive elbow is conserved. 



unity) to: 

P{X - Y = 9) (X ^/[m(ix{gi, - 9)] - mm{gi+i, hi+i - 6)]'^ + [max(6' -I- gi, hi) - mm{6 + gi+i,hi+i]^ . (32) 
This is easily obtained geometrically (the square root conies from Pythagoras' theorem) and the problem results in 



finding the intersect of a Une with the grid. There are two configurations. Working out the cases shows that Eq. (32) 
reduces to a triangular or a top-head truncated triangular distribution. The coefficients entering / can then be taken 
as the potentials weighted by the area intersecting their corresponding region. The soundness of such an approach 
can be appreciated only if cases are checked numerically to compare quantitatively various parametrizations. For 
simplicity, we shall here consider cases where only the dominant potential is considered. An example of such a gap 
equation (/ — 1)A = is defined with: 

Voli V1I2 

^ili K)l3 I , (33) 




where 



tanh(e/(2fcBTc))^^^ ^3^^ 



andTj is integrated on the respective steps, as defined by the integral boundary conditions, i.e., (wf ,a;^,a;j',a;;^,ci;^) 
(0,tJD,WD,a;c,wc,WB). 
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FIG. 13. (solid blue) Three-step potential that models the form of effective electron-electron interaction when polariton- 
mediated with the three parameters ciJb,c,d and, (thick purple) the gap A(^) as computed numerically for this potential. In 
good approximation, the gap is also itself a three-steps function. 



As before, we estimate the critical temperature by the condition that gives vanishing values of the gap on the whole 
interval. The first integral Ii is evaluted in the same approximation as for the usual BCS (or Bogoliubov) potential: 



(35) 



I2 and I3 are also essentially logarithmic in their energy range and we take: 

I2 = ln(a;c/'^D) , and I3 = ln(a;B/wc) ■ (36) 
We solve the gap equation by setting its determinant to zero and solving for Tq, giving 



^bTc = 1.13Swd exp 



Vo 



\ 



1 + 



l+ln 



Vo 



(37) 



We see that the impact of the polariton potential shape on the gap is rather intricate, here as well, with some 
unexpected behaviour, produced both numerically and from this formula. In Fig. |14[ we show the effect of widening the 
central, attractive region of the potential. Naively one would expect this to increase the gap (or critical temperature), 
since A(0) increases linearly with huD. However, Eq. (37) predicts a decrease of the critical temperature, as shown 
in Fig. |15[ The reason why is understood by computing the gap at zero temperature, as seen in Fig. |14[ where the 
repulsive barrier is shown to behave as a trap for the gap function, which results, as the width is increased, in a 
loosening of the gap strength, akin to the quantum mechanical situation of a bound state in a square potential. The 
numerical and linearized models display a good qualitative agreement, as seen in Fig. |15| Let us repeat that we are 
not considering here quantities that can be directly compared, since A(Tc) should be used for this purpose. 

Another example is shown in Fig. [T6j where both the attractive and repulsive parts of the potential are increased 
together. This results in a rapid increase of the gap. The result is also confirmed by numerical simulations, as shown 
in Fig. |17[ and a good agreement is obtained. Unexpectedly, if the repulsive part grows twice as fast as the attractive 
one, not only does this also result in an increase of the gap, but even a faster one. This, too, is confirmed by the 
numerics. Note that in this case, we reach a region where our numerical procedure jumps to other solutions (usually 
of a highly oscillatory character, such as those reported in |55J)j but behaves as expected until then. In this regard, 
the value of having an approximate analytical solution is obvious. 

This analysis is finally applied to the case of the gap equation with the numerical results of the critical temperature 
obtained in both the polariton and exciton cases for the cases of Figs. [6] and [7| The parameters are gathered in 
Table |T] and the results plotted in Fig. [18] Both cases show a strong variation of the critical temperature with 
moderate variations of the condensate density (one order of magnitude). The polariton case is steeper and roughly 
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FIG. 14. Potential (filled blue) and its corresponding gap equation, solved numerically, as the size of the central attractive 
plateau is increased. Somewhat unexpectedly, the gap decreases as a result. The repulsive barriers have the effect of squeezing 
the gap up. 
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FIG. 15. A(0) (joined points, blue) and k^Tc from the three-step Bogoliubov approximation, in the case where the width of 
the central attraction is increased. Unexpectedly, this results in a decrease of the critical temperature. 



linear while the exciton case increases less quickly but begins earlier. In both cases, temperatures are very high, 
and other effects will surely break the mechanism, for instance the loss of the condensate. This shows, however, the 
robustness of the mechanism in the conditions where it should apply. 



V. CONCLUSIONS 



We have studied possible mechanisms of superconductivity in semiconductor heterostructure systems where a Bose- 
Einstein condensate mediating the effective electron-electron interaction leads to enhancement of the coupling, yielding 
very high critical temperatures. We have considered the case of an exciton BEC, consisting of a sandwich of an n- 
doped QW containing the superconducting electrons, in contact with coupled QWs where the BEC of indirect excitons 
is formed, for instance by optical excitation. We have also considered the case of a polariton BEC, where n-doped 
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FIG. 16. Numerical solution of the gap equation when both the attractive and repulsive parts of the potential are increased 
together (here in the same proportions), showing a rapid increase of the gap. 
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FIG. 17. A(0) (joined points, blue) and ksTc from the three-step Bogoliubov approximation. 
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FIG. 18. (colour online) Critical temperature as a function of the condensate density A^o for the case of a polariton condensate 
(solid blue) and of an exciton condensate (dotted purple). 



QWs and undoped QWs hosting excitons are embedded in a microcavity. We have computed in these two cases the 
effective electron-electron interaction U{io) as a function of the exchanged energy huj, showing how the retardation 
effect acquires a peculiar character in the polariton case, namely, yielding a weakly attractive potential at long 
times, followed by a succession of strongly attractive and strongly repulsive windows. The gap equation in this case 
exhibits strong differences as opposed to the case of the Cooper potential (of conventional BCS and exciton-BEC 
mechanism). To understand the physical mechanism leading to large gaps and bypass numerical instabilities, we 
studied in detail the gap equation in the cases of simplified step-wise potentials, and offered an analytical method 
to obtain the critical temperature, which is in qualitative agreement with the numerical results. Our results suggest 
record breaking critical temperatures in these systems. We stress, however, that the path towards achievement of 
exciton mediated superconductivity at high temperatures may be long and full of obstacles. Of the two experimentally 
relevant systems that we have considered, one (coupled QWs) only shows the EEC of excitons at very low temperatures 
(less than IK |38J ) . Consequently, one cannot expect high temperature superconductivity in this system, while at low 
temperatures the superconducting gap may be very large. In microcavities, polariton BEC or polariton lasing have 
indeed been demonstrated at room temperature |17l I18j. However, embedding a high-quality n-doped QW inside 
the cavity, the proper choice of the experimental geometry in order to minimise optical absorption in the doped 
QW, and especially fabrication of quantum contacts for selective injection of carriers in the QW of interest, may 
pose technological difficulties. It is possible that hybrid metal-semiconductor or semi-metal semiconductor systems 
may appear more suitable for the observation of the predicted effects. A conventional superconductor put in contact 
with a semiconductor also seems promising. Collective quantum phenomena in Bose-Fermi mixtures are extremely 
complicated and we foresee breakthroughs in their study in multilayer structures combining a Fermi gas of electrons 
and an exciton BEC. 
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